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Abstract: 

In this article we present a method to implement orthogonal polynomials and many other special 
functions in Computer Algebra systems enabling the user to work with those functions appropriately, 
and in particular to verify different types of identities for those functions. Some of these identities 
like differential equations, power series representations, and hypergeometric representations can even 
dealt with algorithmically, i. e. they can be computed by the Computer Algebra system, rather than 
only verified. 

The types of functions that can be treated by the given technique cover the generalized hypergeometric 
functions, and therefore most of the special functions that can be found in mathematical dictionaries. 

The types of identities for which we present verification algorithms cover differential equations, power 
series representations, identities of the Rodrigues type, hypergeometric representations, and algorithms 
containing symbolic sums. 

The current implementations of special functions in existing Computer Algebra systems do not meet 
these high standards as we shall show in examples. They should be modified, and we show results of 
our implementations. 

1 Introduction 

Many special functions can be looked at from the following point of view: They represent functions 
f(n,x) of one "discrete" variable n G D defined on a set D that has the property that n G D =>• 
n + 1 G D (or n G D =^ n — 1 G D), e. g. D = No,2,H, or (D, and one "continuous" variable 
x £ I where / represents a real interval, either finite / = [a, b], infinite (/ = [a, oo), / = (— oo, a], or 
/ = R), or a subset of the complex plane <D. 

In the given situation we may speak of the family (f n )neD of functions f n {x) := f(n,x). 

In this paper we will deal with special functions and orthogonal polynomials of a real/complex 
variable x. Many of our results can be generalized to special and orthogonal functions of a discrete 
variable x which we will consider in a forthcoming paper. 

Many of those families, especially all families of orthogonal polynomials, have the following 
properties: 
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1. (Derivative rule) 

The functions f n are differentiable with respect to the variable x, and satisfy a derivative rule 
of the form 

q m—1 m—1 

fn( x ) = jrfn{x) = r k(n,x) f n -k(x) Or f' n {x) = ^ r k (jl, x) f n+k (x) , (1) 

OX k=0 k=0 

where the derivative with respect to x is represented by a finite number of lower or higher 
indexed functions of the family, and where r k are rational functions in x. If r m _i(n, x) ^ 
then the number m is called the order of the given derivative rule. We call the two different 
types of derivative rules backward and forward derivative rule, respectively. 

2. (Differential equation) 

The functions f n are m times differentiable (m 6 N) with respect to the variable x, and 
satisfy a homogeneous linear differential equation 

m 

J2Pk(n,x)fi h \x) = 0, (2) 

k=0 

where p k are polynomials in x. If p m (n, k) ^ then the number m is called the order of the 
given differential equation. 

3. (Recurrence equation) 

The functions f n satisfy a homogeneous linear recurrence equation with respect to n 

m 

^2q k (n,x) f n _ k (x) = , (3) 

k=0 

where q k are polynomials in x, and m £ IN". If qo(n,k),q m (n,k) ^ then the number m is 
called the order of the given recurrence equation. 

Some of those families, especially all "classical" families of orthogonal polynomials, have the fol- 
lowing further property: 

4. (Rodrigues representation) 

The functions f n have a representation of the Rodrigues type 

1 d n 

UX) = K^{x)dx-^nix) (4) 

for some functions g depending on x, and h n depending on n and x, and a constant K n 
depending on n. 

From an algebraic point of view these properties read as follows: Let K[x] denote the field of 
rational functions over K where K is one of Q, H, or (D. Then if the coefficients of the occurring 
polynomials and rational functions are elements of K, 

1. the derivative rule states that f n is an element of the linear space over K [x] which is generated 
by {/„, fn-i, /n-(m-i)} or {fn, fn+i, • • • , fn+m-i}, respectively; 
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2. the differential equation states that the m + 1 functions fn (k = 0, . . . , m) are linearly 

!k) 

dependent over if [x]; moreover, by an induction argument, any m + 1 functions f n {k € INo) 
are linearly dependent over if [a;]; 

3. the recurrence equation states that the m + 1 functions f n -k {k = 0, . . . , m) are linearly 
dependent over if [x]; moreover, by an induction argument, any m + 1 functions /„ (n € D), 
are linearly dependent over if [x]. 

One important question when dealing with special functions is the following: Which properties of 
those functions does one have to know to be able to establish various types of identities that those 
functions satisfy? With respect to the implementation of special functions in Computer Algebra 
systems this question reads: Which properties should be implemented for those functions, and in 
which form should this be done such that the user is enabled to verify various types of identities, 
or at least to implement algorithms for this purpose? 



Nikiforov and Uvarov [18] gave a unified introduction to special functions of mathematical 
physics based primarily on the Rodrigues formula and the differential equation. They dealt, how- 
ever, only with second order differential equations, which makes their treatment quite restricted, 
and moreover their development does not have algorithmic applications. 



Truesdell 25] gave a unified approach to special functions based entirely on a special form of 
the derivative rule. His development has some algorithmic content, which, however, is difficult 
or impossible to implement in Computer Algebra. Truesdell's approach — although nice — has the 
further disadvantage that one can obtain only results of a very special form, see |l3j . 

From the algorithmic point of view another approach is better: We will base our treatment of 
special functions on the derivative rule (|l]) in combination with the recurrence equation (g). We will 
show that an implementation of special functions in Computer Algebra systems based on these two 
properties gives a simplification mechanism at hand which, in particular, enables the user to verify 
many kinds of identities for those functions. Some of these identities like differential equations, and 
power series representations can even be dealt with algorithmically, i. e. they can be computed by 
the Computer Algebra system. 

Our treatment is connected with the holonomic system approach due to Zeilberger [|27|]— 1|29|| 
which is based on the valididy of partial differential equations, mixed recurrence equations, and 
difference-differential equations. This connection will be made more precise later. 

The class of functions that can be treated this way contains the Airy functions Ai (x), Bi (x) 
(see e. g. Q, § 10.4), the Bessel functions J n (x),Y n (x), I n (x), and if n (x) (see e. g. [||, Ch. 9-11), 
the Hankel functions Hn (x) and Hn(x) (see e. g. Q, Ch. 9), the Kummer functions M(a,b,x) = 

and U(a,b,x) (see e. g. ||, Ch. 13), the Whittaker functions M n ^ m (x) and W ntTn (x) 

§ 13.4), the associated Legendre functions P%(x) and Q a {x) (see e. g. ||, § 8), all kinds 
of orthogonal polynomials: the Jacobi polynomials P^ a (x), the Gegenbauer polynomials C n a \x), 
the Chebyshev polynomials of the first kind T n (x) and of the second kind U n (x), the Legendre 
polynomials P n (x), the Laguerre polynomials L^\x), and the Hermite polynomials H n (x) (see 
p3| , [24], and Q, § 22), many more special functions, and furthermore sums, products, derivatives, 
antiderivatives, and the composition with rational functions and rational powers of those functions 
(see p|, p7l], pi and E3). 



(see e. g. 



In the case of the classical orthogonal polynomials the properties above can be made much more 
precise (see e. g. p4j, Kapitel IV). Therefore let f n : [a,b] — > K (n e No) denote the family of 
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orthogonal polynomials 

fn{%) = k n x + k n x + . . . 
with respect to the weight function w(x) > 0, i. e. with the property that 

b 

w(x) f n (x) f m (x) dx = (n^m) 



and 

b 

p2 



w(x) f n (x) dx = h n j^0 . 



Then we have the properties: 



1. (Derivative rule) 

The functions f n satisfy a derivative rule of the form 



X f n = (i n f n -l + (^ X " x + a n) fn 



(see e. g. [24j, p. 135, formula (4.8)) where 

a n = nX(0)--X —, /3 n = -- — IKi k\ — X 

Z K n ftji— 1 K n \ I 

and 

!(b — x)(x — a) if a, b are finite 

x — a if b = oo . (5) 

1 if —a, b = oo 

Especially is the order of the derivative rule 2. 

2. (Differential equation) 

The functions f n satisfy the homogeneous linear differential equation with polynomial coeffi- 
cients 

X f n '(x) + K x h f' n {x) + A n f n {x) = 



(see e. g. [24|, p. 133, formula (4.1)) where 

An = ~n \K\ ki - X"^j , 

and X{x) is given by (|J). Especially is the order of the differential equation 2. 

3. (Recurrence equation) 

The functions f n satisfy the recurrence equation 

fn+l(x) = ~C n f n ^(x) + (A n X + B n ) f n (x) (6) 

(see e. g. [24|, p. 126, formula (2.1)) with 

. _ fcft+i R _ fen+i / k n+1 k n \ n — ^"+1 k n -i h n 

Especially is the order of the recurrence equation 2. 



4 



4. (Rodrigues representation) 

The functions f n have a representation of the Rodrigues type 

(see e. g. pi] , p. 129, formula (3.2)), where X(x) is given by (||), i. e. (|j) is valid with 
g{x) = w(x), and h n (x) = w{x) X(x) n . Especially: The order of the polynomial X{x) is < 2. 

Further it turns out that in the case of classical orthogonal polynomials all coefficient functions of 
f n -k are rational also with respect to the variable n, a fact that depends, however, on the special 
normalizations that are used in these cases. 

We mention that no system of orthogonal polynomials besides the classical ones satisfies a 
Rodrigues representation of type ([?]) with a polynomial X (see e. g. p3|] , Kapitel IV, §3). 

We note that using the recurrence equation @, which is valid also for non-classical orthogonal 
polynomials, or any recurrence equation of type (g) of order two (also called three-term recursion), 
recursively, each (backward or forward) derivative rule (Q) is equivalent to a derivative rule 

fn( x ) = k(n,x) f n (x) +l(n,x) f n+1 (x) (8) 

(k, I rational functions with respect to x) of order two. In general, the order of the derivative rule 
can always be assumed to be less than or equal to the order of the recurrence equation. In some 
nice work 



Truesdell presented a treatment of special functions entirely based on the functional 
equation (3). He showed that this difference-differential equation is independent of the differential 
equation (2) and the recurrence equation (||), i. e. it does not imply the existence of one of these. 
In contrast to this work, our main notion is the 

Definition (Admissible family of special functions) We call a family f n of special functions 
admissible if the functions f n satisfy a recurrence equation of type @ and a derivative rule of type 
(|l]). We call the order of the recurrence equation the order of the admissible family /„. □ 

Note that the recurrence equation (g) together with m initial functions f m , f m +i, ■ ■ ■ , fno+m-l 
determine the functions f n (n £ D) uniquely. 

So an admissible family of special functions (with given initial functions) is overdetermined by 
its two defining properties, i. e. the recurrence equation and the derivative rule must be compatible. 
This fact, however, gives our notion a considerable strength: 

Theorem 1 For any admissible family f n of order m the linear space Vf n over K[x] of functions 

generated by the set of shifted derivatives {fn± k \ j, k S No} is at most m-dimensional. On the 

other hand, if the family \ j,k £ No} spans an m-dimensional linear space, then /„ forms an 

admissible family of order m. 

Proof: By the recurrence equation and an induction argument it follows that the linear space V 
spanned by {f n ±k I k £ No} is at most m-dimensional. Using the derivative rule, by a further 
induction it follows that the derivative of any order fjf^ (k € No) is an element of V. Therefore 

v fn = v. 

If on the other hand for a family f n the set of derivatives {f^±k \ 3>k £ No} is m-dimensional, 
then the existence of a recurrence equation and a derivative rule of order m are obvious. □ 
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From the algebraic point of view this is the main reason for the importance of admissible families: 
Any m + 1 distinguished elements of Vf n are linearly dependent, i. e. any arbitrary element of Vf n 
can be represented by a linear combination (with respect to -^[a?]) of any m of the others. This 
is the algebraic background for the fact that so many identities between the members and their 
derivatives of an admissible family exist. 
In particular we have 

Corollary 1 Any admissible family /„ of order m satisfies a simple differential equation of order 
m. □ 

In § H we give an algorithm which, in particular, generates this differential equation of f n . 

With regard to Zeilberger's approach Corollary [l] can be interpreted as follows: Any admissible 
family f n {x) forms a holonomic system with respect to the two variables n, and x, whose defining 
recurrence equation, and the differential equation corresponding to Corollary |l] together with the 
initial conditions 

if } (0), and f k (0) (* = 0,...,ro-l) (9) 



yield the canonical holonomic representation of f n (x) (see [27], Lemma 4.1). 

On the other hand, not all holonomic systems f n (x) form admissible families so that our notion 
is stronger: Let f n (x) := Ai (x) for all n S 2, then obviously f n {x) is the holonomic system 
generated by the equations 

fn( x ) = x fn(x) , fn+l(x) = f n {x) , 

and some initial values, that does not form an admissible family as the derivative f' n is linearly 
independent of {f n \ n E Z} over K[x], see § |j| and thus no derivative rule of the form (|l|) exists. 

A further advantage of our approach is the separation of the variables, i. e. the work with 
ordinary differential equations, and one- variable recurrence equations rather than partial differential 
equations, mixed recurrence equations, and difference-differential equations. So our approach — if 
applicable — seems to be more natural. 

To present an example of an admissible family that cannot be found in mathematical dictio- 
naries, we consider the functions 



tt/2 

2 



2 f 

k n (x) := — / cos (x tan 9 — n 9) 
vr J 



d6 



that Bateman introduced in |4j, see also [14]. He verified that formula (2.7)) 



F n (x) := (-l) n k 2n (x) = {-l) n e- x (L n (2x) - L n ^{2x)) . (10) 

We call F n the family of Bateman functions which turns out to be an admissible family of order 
two. 

Bateman obtained the property ([§], formula (4.1)) 

(n - 1) (F n (x) - F n _i(s)) + (n + 1) (F n {x) - F n+X {xj) = 2xF n (x) 

leading to 

n F n (x) - 2 (n - 1 - x) F n -!(x) + (n - 2) F„_ 2 (x) = (11) 
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which is a recurrence equation of type (^) and order two that determines the Bateman functions 
uniquely using the two initial functions 

F (x) = e~ x and F^x) = -2xe~ x 



which follow from (10). 



Bateman obtained further a difference differential equation (GJ, formula (4.2)) 

{n + l)F n+1 {x)-{n-l)F n ^{x) = 2xF' n {x) , (12) 
which can be brought into the form 

F' n (x) = ^-(in-x)F n {x)-(n-l)F n ^ 1 (x)) (13) 

using (|TT[) . This is a derivative rule of the form (|l|) and order two. Therefore F n {x) form an 
admissible family of order two. 

We note that the functions F n satisfy the differential equation 

x F' n '{x) + {2n - x) F n {x) = 0, (14) 

(see formula (5.1)), and the Rodrigues type representation 

(see g|, formula (31)). 

2 Properties of admissible families 

Theorem 2 Let f n form an admissible family of order m. Then 

(a) (Shift) f n ±k (k £ IN) forms an admissible family of order m; 

(b) (Derivative) f' n forms an admissible family of order < m; 

(c) (Composition) f n o r forms an admissible family of order < m, if r is a rational function, 
and of order < mq, if r(x) = x v l q (p, q € IN). 

If furthermore g n forms an admissible family of order < I, then moreover 

(d) (Sum) /„ + g n forms an admissible family of order < m + I; 

(e) (Product) f n g n forms an admissible family of order < ml. 

Proof: (a): This is an obvious consequence of Theorem |]. 

(b): Let g n ■= fn- We start with the recurrence equation for /„ and take derivative to get 

m m 

l'k( n > x ) fn-k{x) + J2lk(n,x)f n _ k (x) = . (16) 

k=0 k=0 
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From Theorem |l], we know that each of the functions f n —j (j = 0, . . . , m) can be represented as a 
linear combination of the functions f' n __ k (k = 0, . . . ,m — 1) over if[x] 3 which generates a recurrence 
equation for g n . Similarly a derivative rule for g n is obtained. 

(c) : For the composition h n := f n ° r with a rational function r, the recurrence equation is obtained 
by substitution, and the derivative rule is a result of the chain rule. If r(x) = x l l q , then, by |l5| , 
Lemma 1, the family {hn \ j G INo} is spanned by the mq functions x r l q 0{x 1 / q ) (j = 1, . . . , m — l, 
r = 0, ... ,q — 1), and since {f^±k I 3i^ e ^o} has dimension m, the linear space spanned by 

{h^lj. I 3i^ ^ ^o} nas dimension < mq, implying the result. If finally r(x) = x p ^ q , then a 
combination gives the result. 

(d) : By a simple algebraic argument, we see that f n -k + g n ~k (k G Z) span the linear space 
V := Vf n+9n = Vf n + V^ n of dimension < m + Z over ET[x]. Therefore /„ + g n satisfies a recurrence 
equation of order < m + I. If we add the derivative rules for f n and g n , we see that f n + g' n G V, 
and thus can be represented in the desired way. 

(e) : By a similar algebraic argument (see e. g. [^], Theorem 2.3) we see that f n -k ■ 9n~k {k G 2) 
span a linear space of dimension < ml over if [a;], hence f n g n satisfies a recurrence equation of 
order <ml. By the product rule, and the derivative rules for /„ and g n we see that the derivative 
of fn9n is represented by products of the form f n ~k9n-j (k,j G 2), and as those span the linear 



space V (see e. g. |15[, Theorem 3 (d)), we are done. □ 



As an application we again may state that the Bateman functions form an admissible family: Using 
the theorem, this follows immediately from representation (|io|). 

Next we study algorithmic versions of the theorem. The following algorithm generates a repre- 
sentation of the members f n ±k (k = 0, . . . ,m — 1) of an admissible family in terms of the derivatives 
fn±j (j = 0, . . . , m — 1). By Theorem [l] we know that such a representation exists. Without loss of 
generality, we assume that the admissible family is given by a backward derivative rule. In case of 
a forward derivative rule, a similar algorithm is valid. 

Algorithm 1 Let f n be an admissible family of order m, given by a backward derivative rule 

m—l 

fn( x ) = r k( n ,x) fn-k(x) • 
fc=0 

Then the following algorithm generates a list of backward rules (fc = 0, . . . , m — 1) 

m—l 

f n . k (x)= X^JCn,*)/^) (17) 

3=0 



(Rj rational with respect to x) for f n -k (k = 0, . . . , m — 1) in terms of the derivatives f n _j (j 
0,...,m-l): 

(a) Shift the derivative rule m — l times to obtain the set of m equations 

m—l 

fli-j( x ) = r k( n ~ 3, x) / n _j_ fc (x) (j = 0, . . . , m - 1) . 



k=0 
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(b) Utilize the recurrence equation to express all expressions on the right hand sides of these 
equations in terms of f n -k (k = 0, . . . , m — 1) leading to 

m— 1 

f' n -j{x) = 22, r k( n > x ) fn-k(x) (j ' = 0, . . . , m — 1 , rj, rational with respect to x) . 

(c) Solve this linear equations system for the variables (k = 0, . . . , m — 1) to obtain the 
representations (|l7"|) searched for. □ 

The proof of the algorithm is obvious. It is also clear how the method can be adapted to obtain 
forward rules in terms of the derivatives. As an example, the algorithm generates the following 
representations for the Bateman functions 

F n {x) = 1 ~ " + X F' n {x) + W " 1 F' n _ x {x) , 
In — 1 — x ln—1 — x 



and 

1 "I - 72 x j . . 1 + n , 



n+1 



1 + 2n - x l + 2n — x 

in terms of their derivatives. 

We note that by means of Algorithm | and the results of ||] (see also P- 342, and HU), 
we are able to state algorithmic versions of the statements of Theorem 

Algorithm 2 The following algorithms lead to the derivative rules and recurrence equations of 
the admissible families presented in Theorem ^ 

(a) (Shift) Direct use of derivative rule and recurrence equation lead to the derivative rule and 
the recurrence equation for f n ±\; a recursive application gives the results for f n ±k (k £ IN"). 



(b) (Derivative) By Algorithm [l| we may replace all occurrences of f n -k (k = 0, . . . , m) in (|l~6|), 
resulting in the recurrence equation for f n \ similarly the derivative rule is obtained. 

(c) (Composition) If r is a rational function, then an application of the chain rule leads to the 
derivative rule and the recurrence equation of f n o r; an approach similar to the algorithmic 
version of Theorem 2 in |l5| yields the derivative rule and the recurrence equation of f n o x 1 / 9 
by an elimination of the expressions x r / q fn \x 1 / q ) (r = 1, . . . , q — 1, j = 1, . . . , m — 1). 



(d) (Sum) Applying a discrete version of Theorem 3 (c) in |15| to /„ + g n (see also |27|| , p. 342, 
and |^T|, Maple function rec+rec) results in the recurrence equation, and a similar approach 
gives the derivative rule. 

(e) (Product) Applying a discrete version of Theorem 3 (d) in IS] to f n g n (see also [27], p. 342, 
and plj, Maple function rec*rec) yields the recurrence equation, and a similar approach 
gives the derivative rule. □ 

A Mathematica implementation of the given algorithms generate e. g. for the derivative F' n (x) of 
the Bateman function F n (x) the derivative rule 

F "( x ) = o — -2 ( - 1) CiW + (1 " n + x) F n (x 

x — 2nx + x zx - 
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and the recurrence equation 

K+i(x) = n , w. 1 9 - v ((n-l)(x-2n-l)F;_ 1 (x) + 2 (l-2n 2 +3nx-x 2 )F^(x) 
{l+n)(l— 2n+x) \ 

and for the product A n {x) := -F 2 (x) the derivative rule 

At I \ (1 ~~ n ) ( n ~~ / \ 

= ^ 77^ 7^^-21 

2 n x (1 — n + x) 

+ 2(n-l) + 

n x 

(3n - 3n 2 - 4x + 8nx - 4x 2 ) 
2x (1 — n + x) 

and the recurrence equation 

_ 1 f (n-2) 2 (n-l)(s-n) 
^ n+UXj " (1+n) 2 \ n(l-n + x) A »" 2 

(re - 1) (3 re - 3n 2 - 4x + 8nx - 4x 2 ) „ 

+ K-i 

re 

(x — n) (—3 n + 3 n 2 + 4 x — 8 n x + 4 x 2 ) \ 

+ : l A n 

1 — n + x I 

are derived. 



3 Derivative rules of special functions 



Many Computer Algebra systems like Axiom |J, Macsyma |16| , Maple [pTj , Mathematica p6| , 
or Reduce ^ support the work with special functions. On the other hand, there are so many 
identities for special functions that it is a nontrivial task to decide which properties should be used 
by the system (and in which way) for the work with those functions. 

Since all Computer Algebra systems support derivatives, as a first question it is natural to ask 
how the current implementations of Computer Algebra systems handle the derivatives of special 
functions. Here are some examples: Mathematica (Version 2.2) gives 

In[l] := D[BesselI[n,x] ,x] 

Bessell[-1 + n, x] + Bessell [1 + n, x] 



Out[l]= 

2 

In[2] := D[LaguerreL[n,a,x] ,x] 

Out [2]= -LaguerreL [-1 + n, 1 + a, x] 

We note that in Mathematica the derivatives of all special functions symbolically are imple- 
mented. On the other hand, we notice that, given the function I n (x), Mathematical derivative 
introduces two new functions: l n -\ (x), and I n +i (x). Given the Laguerre polynomial Ln(x), the 
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derivative produced introduces a new function where both n, and a are altered. The representation 
used is optimal for numerical purposes, but is not a representation according to our classification. 
With Maple (Version V.2) we get 

> diff (Bessell(n,x) ,x) ; 

n Bessell(n, x) 

BesselKn + 1, x) + 

x 



> dif f (L(n, a,x) ,x) ; 

d 

L(n, a, x) 

dx 

Thus Maple's derivative for the Bessel function I n (x) introduces only one new function / n +i(x), 
and is of type (Q), whereas (even if orthopoly is loaded) no symbolic derivative of the Laguerre 
polynomial L n °'\x) is implemented. 

Obviously there is no unique way to declare the derivative of a special function. However, we 
note that if we declare the derivative of a special function by a derivative rule of type (|l[) of order 
m then we can be sure that the derivative of the special function f n (x) introduces at most m new 
functions, namely f n -k{x) (k = 1, . . . ,m). Moreover, if the family of special functions depends on 
several parameters, then the given representation of the derivative does not use any functions with 
other parameters changed. 

Here we give a list of the backward derivative rules of the form ([[]) for the families of special 
functions that we introduced in § |l| which all turn out to be of order two (see e. g. Q, (9.1.27) 
(Bessel and Hankel functions), (9.2.26) (Bessel functions), (13.4.11), (13.4.26) (Kummer functions), 
(13.4.29H13.4.33) (Whittaker functions), (8.5.4) (associated Legendre functions), and § 22.8 (or- 
thogonal polynomials)): 

- Jn (x) 

x 

-Y n (x) 

x 

- In {x) 

X 

n 

K n 

x 

d „m, s „m , . n 

dx " ' .r 

. [■■») ri 

dx " ' ' '" ~ " ' x 

—M(a,b,x) = -((b - a)M(a -l,b,x) - (b- a - x)M(a,b,x)) 
ox x V / 

— U (a,b,x) = — ( — U(a — 1, 6, x) + (a — b + x) U(a, b,x)) , 
ox ' x V ' / 



J'nix) 


= J n -i (x) - 


Y>(x) 


= Y n ^{x)- 


I'n ( x ) 


= In-! (x) - 


K(x) 


= -K n -i (x 


2#>(z) 


= n&ix)- 


H n 2 \x) 


= H^ix)- 
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M' nm {x) = ^((l + 2m-2n)M n _ lim (x) + (2n-x)M n>m (x 
W^ m (x) = — ((1 - 4m 2 - 4n + 4n 2 ) W n - 1>m (x) + (4n - 2x) W n , m (x)) , 



<9x 



-^ P t P X*) = {2n+a+ 1 m _ x 2 ) (2(n+a){n+P)P^f (x) + n(«-/3- (2n+a+)9)x)P^ (x)) , 

CW(x) = T ^ f ((n + 2a-l)C't ) 1 (x)-nxC( Q )(x)) , 

(x) = - — (reT n _i (x) -nxT n (xfj , 

U' n {x) = - — ^ 2 ((n + 1) {7 n _i (x) - n x Z7 TO (x)) , 



9a 
<9x 



1 J*" 2 in-i (x) - n x P n (x)) , 



— 4«)(x) = i(-(n + a)LS2 l(a . )+nL W (a .)) , (18) 
iln(x) = 2nH n -i(x) . 

4 Recurrence equations of special functions 

Whenever in any expression subexpressions of the form rfc / n _fc (r^ rational, k G Z) occur, in an 
admissible family of order m with the recursive use of the recurrence equation we may replace so 
many occurrences of those expressions f n -k that finally only m successive terms of the same type 
remain. 

This allows for example to eliminate the number of occurrences in any linear combination (over 
-fT[x]) of derivatives of f n to m, a fact with which we will deal in more detail in § ||. 

We show how Mathematica and Maple work with regard to this question. Whereas Math- 
ematica does not have any built-in capabilities to simplify the following linear combinations of 
Bessel and Laguerre functions, 

In [3] := Bessell [n+l,x]+2*n/x*BesselI[n,x] -Bessell [n-l,x] 

2 n Bessell [n, x] 

Out [3]= -Bessell [-1 + n, x] + + Bessell [1 + n, x] 

x 
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In [4] := Simplify [7J 

2 n Bessell [n, x] 

Out [4]= -Bessell [-1 + n, x] + + Bessell [1 + n, x] 

x 

In [5] : = LaguerreL [n+1 , a , x] - (2*n+a+l-x) *LaguerreL [n , a, x] + (n+a) *LaguerreL [n-1 , a, x] 
Out [5] = (a + n) LaguerreL [-1 + n, a, x] - 

> (l + a + 2n-x) LaguerreL [n, a, x] + LaguerreL [1 + n, a, x] 
In [6] := Simplify ['/.] 

Out [6] = (a + n) LaguerreL [-1 + n, a, x] - 

> (l + a + 2n-x) LaguerreL [n, a, x] + LaguerreL [1 + n, a, x] 

with Maple we get 

> BesselI(n+l,x)+2*n/x*BesselI(n,x)-BesselI(n-l,x) ; 

n Bessell (n, x) 

BesselKn + 1, x) + 2 BesselK- 1 + n, x) 

x 

> simplifyC') ; 



> L(n+l,a,x)-(2*n+a+l-x)*L(n,a,x)+(n+a)*L(n-l,a,x) ; 

L(n + 1, a, x) - (2 n + a + 1 - x) L(n, a, x) + (n + a) L(n - 1, a, x) 

> simplifyC) ; 

L(n + 1, a, x) - 2 L(n, a, x) n - L(n, a, x) a - L(n, a, x) + L(n, a, x) x 
+ L(n - 1, a, x) n + L(n - 1, a, x) a 

i. e. Maple's simplify command supports simplification with the aid of the recurrence equations 
for the Bessel functions. On the other hand, for the orthogonal polynomials (even if orthopoly is 
loaded) no simplifications occur. 

In the rest of this section we give a list of the recurrence equations of the given type for the 
families of special functions that we consider which all turn out to be of order two (see e. g. Q, 
(9.1.27), (9.2.26), (13.4.1), (13.4.15), (13.4.29), (13.4.31), (8.5.3), and § 22.7). We list them in the 
form explicitly solved for F n+ \ as this is the usual form found in mathematical dictionaries. 

In 

J n+ i (x) = -J n -\ (x) H J n (x) , 

x 

2n 

Y n+1 (x) = -y n _i (x) + —Y n (x) , 

x 
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rr(2) 
11 n+l 

M(a + 1 
U{a + l,b. 



r n+l {x 

(x 
(x 

(x 
b, x 



n+l 



», X 



M n+1 , m (x 
W n+1>m (x 

Pa+l (X 
Qa + 1 (X 

T n+ l (x 
U n +l {X 

Pn+1 (X 

l { :Ux 



2n 



= I n -i (x) I n (x) , 

2n 

= K n -i (x) H K n (x), 



x 



2n 



-H n 1 l 1 (x) + ^H^ (x), 
+ (x), 
a)M(a - 1,6, a;) + (2a - 6 + x) M(a, 6, x)) , 
- a^ + a-b) ^ -hb,x) + (b-2a-x)U (a, 6, x)) , 

1 + 2m + 2n ( (1 + 2m " 2n ) M ™~1,™ ( x ) + ( 4n ~ 2x ) M ™>™ ( x )) . 
- ((-1 + 4m 2 + 4n - 4n 2 ) W n _ ljm (x) - (8n - 4x) W n , m (x)) , 

_1_ ( - (a + 6) (x) + (2a + 1) x P^ (x)) , 
_i_ ( - (a + 6) Qt i 0*0 + (2a + 1) x Q*_i (*)) - 

2 (n+l) (n + a + ^ + l) (2n +a+/3 ) (-2(n + a)(n + /3)(2n + a + /? + 2) i£f> (x) 
+ ((2n+a+/?+l)(a 2 -/? 2 ) + (2n+a+/3) 3 x)pi a ' /3) (x)) , 

_1_ (_ ( n + 2a - 1) c£\ (x) + 2(n + a) x C,W (x)) , 

-T n _i(x) + 2xT n (x) , 
-C/ n _i (x) + 2xf7 n (x) , 

— i—(-nP„_i(x) + (2n+l)xP n (x)) , 
n + 1 v / 

^ (-(n + a) L^(x) + (2n + a + 1 - x) 4°)(x)) , 

-2nH n ^i{x) + 2xH n (x) . 



H n+ i{x) 

Note that (a)k (which is used in the recurrence equation for the Jacobi polynomials P^^\x)) 
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k 

denotes the Pochhammer symbol (or shifted factorial) denned by (o)fc := II ( a +J — !)• 

i=i 

We note further that for functions with several "discrete" variables it may happen that for each 
of them there exists a recurrence equation. As an example we consider the Laguerre polynomials 
for which we have (@ (22.7.29), in combination with (22.7.30)) 

L^{x) = ^(-{n + a)L^- 1 \x) + {a + x)L^{x)) . (19) 

In § [7| we will demonstrate that generalized hypergeometric functions satisfy recurrence equations 
with respect to all their parameters. 

To be safely enabled that the algorithms of § |9|-§ [ll] apply, all of those recurrence equations 
should be implemented and applied recursively for simplification purposes. 

5 Embedding of one-variable functions into admissible families 

In this section we consider first, how the elementary transcendental functions are covered by the 
given approach. 

Consider the exponential function f(x) = e x . This function can be embedded into the admissible 
family f n , defined by the properties 

fn( x ) = fn(x) , fn+l(x) = f n {x) and f (x) = e x , 

i. e. the family of iterated derivatives of e x . 

Obviously this is a representation of an admissible family of order one. 

Moreover in the given case it turns out that f n (x) = e x = fo(x) for all u 6 2, so there is no 
actual need to give the functions numbers, and therefore we (obviously) keep the usual notation. 

Similarly the functions sin x and cos x are embedded into the admissible family f n of order two 
given by the properties 

fn(x)=fn-l(x), f n +i(x) = -f n -i(x) , and / (x)=cosx, /i(x) = sinx . 

Again, the family of functions f n is finite, and our numbering is unnecessary: 

if n = 4m (m £ 2) 
if n = 4m + 1 (m £ 2) 
if n = 4m + 2 (m £ Z) 
if n = 4m + 3 (m £ Z) 

Essentially there are only the two functions cosx, and sinx involved. Note, however, that both 
functions are needed as no simple first order differential equation for sinx or cosx exists. 

Other nontrivial examples of essentially finite admissible families of special functions are formed 
by the Airy functions. Let Ai n (x) = Ai^ n ^ (x), i. e. 

Ai' n (x) = Ai n+ i (x) . 



fn(x) 



cos X 
sinx 

— cos x 

— sinx 
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By the differential equation for the Airy functions (see e. g. M, (10.4)) we have Ai"(x)— x Ai(x) = 0, 
so that from Leibniz's rule it follows that 



Ai n+1 (x) = Ai(" +1 )(x) = (Ai"(x)) 



(n-l) 



n-l 



= (x Ai (x)j = X! ( fc ) x(fe) ( Ai ' 

k=0 

= x Ai^" 1 ) (x) + (n - 1) Ai( n " 2 ) (x) = x Ai n _! (x) + (n - 1) Ai n _ 2 (x) , 
and therefore Ai (x) is embedded into the admissible family Ai n of order three given by 

Ai' n (x) = Ai n+ i (x) , Ai n+ i (x) = x Ai n _i (x) + (n-l) Ai n _ 2 (x) , (20) 
and we have the initial functions 

Aio (x) = Ai (x) , Aii (%) = Ai' (x) and Ai 2 (x) = x Ai (x) . 
Similarly Bi (x) is embedded into the admissible family of order three given by 

Bi' n (x) = Bi n+ i (x) , Bi n+ i (x) = x Bi n _i (x) + (n - 1) Bi n _ 2 (x) , (21) 
and the initial functions 

Bio (x) = Bi (x) , Bii (x) = Bi' (x) and Bi 2 (x) = x Bi (x) . 

Our indexed families turn out to be representable by 

Ai n (x) = p n (x) Ai (x) + q n (x) Ai' (x) and Bi n (x) = p n (x) Bi (x) + q n (x) Bi' (x) , 

with polynomials p n and q n in x. This shows, however, that to deal with the Airy functions 
algorithmically as is suggested in this paper, besides the functions Ai (x) and Bi (x) the two 
independent functions Ai' (x) and Bi' (x) are needed, but none else. Let's look, how Computer 
Algebra systems work with the Airy functions. 
Maple handles them as follows: 



> Ai(x) ; 

> diff (Ai(x) ,x) ; 



Ai(x) 



1/2 3/2 
2 BesselK(l/3, 2/3 x ) 



1/4 



1/4 
x Pi 

/ 3/2 \ 

1/2 5/4 I 3/2 BesselK(l/3, 2/3 x )| 

2x1- BesselK(4/3, 2/3 x ) + 1/2 1 

I 3/2 | 

\ x / 

+ 1/3 

Pi 
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> simplify (diff (Ai(x) ,x$2)-x*Ai(x) ) ; 

1/2 3/2 1/2 3/2 3/2 

1/48 (-3 2 BesselK(l/3, 2/3 x ) - 8 2 BesselK(-2/3, 2/3 x ) x 

1/2 3 3/2 9/4 / 5/4 

+ 16 2 x BesselK(l/3, 2/3 x ) - 48 x Ai(x) Pi) / (x Pi) 

/ 

> diff (Bi(x) ,x) ; 

d 

Bi(x) 

dx 

> diff (Bi(x) ,x$2) ; 

2 

d 

Bi(x) 

2 

dx 

So the derivative of Ai (x) is represented by Bessel functions, whereas the function Ai (x) itself is 
not, and therefore the expression diff (Ai(x) ,x$2) -x*Ai (x) is not simplified. On the other hand 
the derivative of Bi (x) is not a valid Maple function. With Mathematica we get 



In [7] : = 


D [AiryAi [x] ,x] 


Out [7] = 


AiryAiPrime [x] 


In [8] : = 


D [AiryAiPrime [x] ,x] 


Out [8] = 


x AiryAi [x] 


In [9] : = 


D [AiryAi [x] , {x , 2}] -x*AiryAi [x] 


Out [9] = 





In[10] : 


= D[AiryBi[x] ,x] 


Out [10] 


= AiryBiPrime [x] 


In[ll] : 


= D [AiryBiPrime [x] ,x] 


Out [11] 


= x AiryBi [x] 


In[12] : 


= D [AiryBi [x] , {x , 2>] -x*AiryBi [x] 


Out [12] 


= 



Thus we see that in this situation Mathematica does exactly what we suggest: It works with the 
independent functions Ai (x), Ai' (x), Bi (x), Bi' (x), and the derivative rules ( p0| ) and (|2l|). 
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As a further example of an admissible family we consider the iterated integrals 

oo 

erfc„, (x) = J erfc„_i (t) dt 

x 

of the (complementary) error function erfc (x) = 1 — erf (x) = erfco (x) (see e. g. Q, (7.2)) that 
form the admissible family with 

erfc n (x) = - erfc n _i (x) , erfc n+ i (x) = — — — erfc„_i (x) — erfc n (x) , 

2[n + lj n + 1 

and the initial functions 

erfco (x) = erfc (x) , erfci (x) = ■= i^pHx erfc (x) — e~ x ) 

^/7^ V / 



/7T 

(one may also use the initial value function erfc_i (x) = -^=e~ x2 ). In particular, erfcx is embedded 
into an admissible family. 

Maple deals with these functions as suggested: 

> dif f (erf c (n,x) ,x) ; 

- erfc(n - 1, x) 

> simplif y(diff (erf c(n,x) ,x$2)+2*x*dif f (erfc(n,x) ,x)-2*n*erfc(n,x)) ; 



As a final example, we mention another family of iterated integrals, the Abramowitz functions 

oo 

A n (x) := J t n e- t2 - x/t dt 
o 

(see jjj], and @, (27.5)) which form an admissible family with derivative rule 

/ oo \ oo oo 

A ' n (x) = ^[[t n e~ t2 ~ x/t dt )=J-^{ tn e~ t2 - x/t ) dt = -[ t n ~ l e"* 2 "^ dt = -A n ^(x) 



of order one (see (27.5.2)), and recurrence formula 

ft X 

A n+ i{x) = - A n ^i{x) + - A n -2{x) 

of order three (§], (27.5.3)). 

Again, embedded into an admissible family, especially the function Aq{x) = J e~ l ~ x ' t dt is 

o 

covered by our approach. 
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6 Embedding the inhomogeneous case 

Some families of functions are characterized by inhomogeneous differential rules and recurrence 
equations. Examples for this situation are the exponential integrals given by 

oo 

f e xt 

En(x)= J —dt 

1 

(see e. g. [^], (5.1)), and the Struve functions H n (x) and L n (x) (see e. g. [||], Chapter 5), for which 
we have the inhomogeneous properties 

E' n (x) = -E n _i (x) , E n+ i (x) = - E n (x) , 

n n 

(g, (5.1.14) and (5.1.26)), 

H n _i(x) - H n+l( x) = 2H1W - 2n 0F ; ( " n + 3/2) , (22) 

2n T n 
H n _i(x) + H n+ i(x) = — H n (x) + 



x v ' 2«0FT(n + 3/2) 



g, (12.1.9)-(12.1.10)), and 



-L n -x{x) + L n+1 (x) = 2L' n (x) - , (23) 

2 n yJ-K T{n + 3/2) 

2n ^ n 
L n -l(x) - L n+ i(x) = — L n (x) + 



x v ' 2™^Fr(n + 3/2) 

(g, (12.2.4)-(12.2.5)), respectively. Eliminating the inhomogeneous parts (using T(3/2 + n) = 
(1/2 + n) r(l/2 + n)), these examples are made into admissible families with the derivative rules 

E' n (x) = -E n _ x (x) , 



n' n (x) = H n _i(x) H„(x) , (24) 



n 



V n (x) = L n _i(x) L n (x) , (25) 

x 



and the recurrence equations 



E n+ i(x) = ~(xE n -x(x) + (n-l-x)E n (x] 



n 



I / x^ + 2n + An? 

Hn+i(x) = — — (xH n _ 2 (x) + (l-4n)H n _ 1 (x) + H n (x) 

zn + 1 v x 

j / ^2 — 2^ — 4ti^ 

L n+ i(x) = - — — ( - xL n _ 2 (x) - (1 - 4n) L n _i(x) H L n (x) 

An + 1 v x 
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so that the exponential integrals form an admissible family of order two, and the Struve functions 
H n (x) and L n (x) form admissible families of order three. Note that the above derivative rules 
(24)— (25) are not listed in although they are much simpler than the inhomogeneous relations 

After bringing the inhomogeneous rules into the desired form, those families are recognized as 
admissible families, and our method can be applied. 



7 Functions of the hyper geometric type as admissible families 

All functions introduced in this paper are special cases of functions of the hypergeometric type (see 
U ) . In this section we will show that the generalized hypergeometric function p F q defined by 



/ ai a 2 ■■■ a 
ptq h b 2 ■■■ b, 



g 



"h ~h>M>>'<!»)k---(b q ) k ki x ' (26) 



and thus by Theorem || (c) all functions of the hypergeometric type, form admissible families. 
Therefore we first deduce a derivative rule of order two for p F q . 

Let us choose any of the numerator parameters n := a k {k = 1, . . . ,p) of p F q as parameter n. 
Further we use the abbreviations 



From the relation 



n a>2 ■ ■ ■ a p 
h b 2 ■■■ b q 



(n + 1)^ _ n + k 



^A k (n)x 

k=0 



k 



(n) k n 



it follows that 

nA k {n + l) = (n + k)A k (n) . 
Using the differential operator 6f(x) = x f'(x), we get by summation 



CO 



nF n+1 (x) = nJ2M n + 1 )x k = (n + k)J2 A k(n)x k 

k=0 k=0 

oo 

= nF n (x) + J2 kA k( n )x k = nF n (x) + 6F n (x) , 

k=0 

and therefore we are led to the derivative rule 

9F n {x) = n (F n+1 {x) - F n (xj) , or F' n {x) = ~ (F n+1 (x) - F n (x)) . (27) 

Hence we have established that for any of the numerator parameters n := a k (k = 1, . . . ,p) of p F q 
such a simple (forward) derivative rule is valid. 

We note that by similar means for each of the denominator parameters n := b k (k = 1, . . . , q) 
of p F q the simple (backward) derivative rule 

OF n (x) = {n-l)(F n „ x (x)-F n {x)) , or F' n {x) = — - (f„_i(x) - F n {xj) (28) 
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is derived. 

Next, we note that F n satisfies the well-known hypergeometric differential equation 

9(9 + 61 - 1) • • • (9 + b q - l)F n {x) = x(9 + ai)(6 + o 2 ) • • • (9 + a p )F n (x) . (29) 

Replacing all occurrences of 9 in ( f29|) recursively by the derivative rule (|27|) or (|28|), a recurrence 
equation for F n is obtained that turns out to have the same order as the differential equation (f29|). 
i. e. max{p, q + 1}. 

We summarize the above results in the following 



Theorem 3 The generalized hypergeometric function p F q 
derivative rules 



ax a 2 
h b 2 



x satisfies the 



6F n (x) = n(F n+1 {x)-F n {x)" / 
for any of its numerator parameters n := (k = 1, . . . ,p), and 

9F n (x) = (n - 1) (F„_i(x) - F n (x)) 

for any of its denominator parameters n := b/~ (k = 1, . . . and recursive substitution of all 
occurrences of 9 in the hypergeometric differential equation 



9(9 + h - 1) • • • {9 + b g - l)F n (x) = x(9 + ai )(9 + a 2 ) 



+ a p )F n (x) 



generates a recurrence equation of the type (||) of order max{p, q+ 1} with respect to the parameter 
chosen. This recurrence equation has coefficients that are rational with respect to x, and n. In 
particular, p F q forms an admissible family of order max{p, 9+1} with respect to all of its parameters 
a k ,b k . □ 

We note that if some of the parameters of p F q are specified, there may exist a lower order differential 
equation, and thus the order of the admissible family may be lower than the theorem states. We 
note further that this theorem is the main reason for the fact that so many special functions 
form admissible families: Most of them can be represented in terms of generalized hypergeometric 
functions. 



8 Algorithmic generation of differential equations 

In this section we show that the algorithm to generate the uniquely determined differential equation 
of type @ of lowest order valid for / which was developed in Q (see also (l5|]), does apply if / is 
constructed from functions that are embedded into admissible families. 

Algorithm 3 (Find a simple differential equation) Let / be a function given by an expres- 
sion that is built from the functions expx, lnx, sinx, cosx, arcsinx, arctanx, and any other func- 
tions that are embedded into admissible families, with the aid of the following procedures: differ- 
entiation, antidifferentiation, addition, multiplication, and the composition with rational functions 
and rational powers. 

Then the following procedure generates a simple differential equation valid for /: 
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(a) Find out whether there exists a simple differential equation for / of order N := 1. Therefore 
differentiate /, and solve the linear equation 

f'(x) + A f(x) = 

for Aq; i. e. set Aq := Is rational in x, then you are done after multiplication with 

its denominator. 

(b) Increase the order N of the differential equation searched for by one. Expand the expression 

f (N) {x) + A^f^ix) + ■■■ + A f(x) , 

apply the recurrence formulas of any admissible family F n of order m involved recursively to 
minimize the occurrences of F n _^ to at most m successive /c-values, and check, if the remaining 
summands contain exactly N rationally independent expressions considering the numbers 
Aq, Ax, . . . , An_i as constants. Just in that case there exists a solution as follows: Sort with 
respect to the rationally independent terms and create a system of linear equations by setting 
their coefficients to zero. Solve this system for the numbers Aq, A\, . . . , Apf—i- Those are 
rational functions in x, and if there is a solution, this solution is unique. After multiplication 
by the common denominator of Ao, A\, . . . , An-i you get the differential equation searched 
for. Finally cancel common factors of the polynomial coefficients. 

(c) If part (b) was not successful, repeat step (b). 

Proof: Theorem 3 of |15| (compare |22|| ) shows that for / a differential equation of type (^) exists. 
We assume that differentiation is done by recursive descent through the expression tree, and an 
application of the chain, product and quotient rules on the corresponding subexpressions. It is clear 
that the algorithm works for members of admissible families, compare Theorem [j] and Corollary [l]. 
Similarly the algorithm obviously works for derivatives and antiderivatives of admissible families. 
Further it is easily seen that the derivatives of sums, products, and the composition with rational 
functions and rational powers form either sums, or sums of products all of which by a recursive use 
of the recurrence equations involved are represented by sums of fixed lengths, compare Theorem^. 
Thus after a finite number of steps, part (b) of the algorithm will succeed (sharp a priory bounds 
for the resulting orders are given in [jl^]). □ 

We note that from the implementational point of view the crucial step of the algorithm is the 
decision of the rational independency in part (b). If this decision can be handled properly, then 
the proof given in Q shows that the algorithm generates the simple differential equation of lowest 
order valid for /. 

In our implementations, for testing whether some terms are rationally dependent, we divide 
each one by any other and test whether the quotient is a rational function in x or not. This is an 
easy and fast approach which never leads to wrong results, but may miss a simpler solution, which 
in practice, rarely happens. 

Typically this happens, however, for orthogonal polynomials with prescribed n, for which a first 
order differential equation exists. In this case, the recurrence equation hides these rational depen- 
dencies, and in some sense (s. ||, § 7) here it is even advantageous that the rational dependency is 
not realized. 

Another example where our implementations yield a differential equation which is not of lowest 
order is given by 



22 



In [13]:= SimpleDE[Sin[2 x] -2 Sin[x] Cos[x],x] 



Out [13]= 4 F[x] + F"[x] ==0 

This happens because the functions sin (2x) and 2 sin x cos x algebraically cannot be verified to be 
rationally dependent even though they are identical. 

We note that, for elementary functions, we could use the Risch normalization procedure [20] to 
generate the rationally independent terms, but this does not work for special functions. 

Further we note that in case of expressions of high complexity, the use of |jl5[] , Algorithm 2, 
typically is faster. This algorithm, however, in general leads to a differential equation of higher 
order than Algorithm ||[ 

As a first application of Algorithm || we consider the Airy functions Ai n , again, for which the 
Mathematica implementation of our algorithm yields 

In [14] := SimpleDE[AiryAi [n,x] ,x] 

(3) 

0ut[14] = (-1 - n) F[x] - x F'[x] + F [x] == 
i. e. the differential equation 

Ai^ (x) - x Ai' n (x) - (n + 1) Ai n (x) = . (30) 
Similarly, we get for the square of the Airy function 
In [15]:= SimpleDE[AiryAi [x] ~2,x] 

(3) 

Out [15]= -2 F[x] - 4 x F'[x] + F [x] == 

The next calculation confirms the differential equation for the Bateman functions F n (|i~4| ) 

In [16] := SimpleDE [Bateman [n,x] ,x] 

Out [16]= (2 n - x) F[x] + x F" [x] ==0 

Other examples are given with the aid of the iterated integrals of the complementary error function, 
and the Abramowitz functions: 

In[17]:= SimpleDE[Erf c [n,x] ,x] 

Out [17]= -2 n F[x] + 2 x F'[x] + F"[x] ==0 

(see § (7.2.2)) and 

In [18] := SimpleDE [Exp [a x]*Erfc[n,x] ,x] 
2 

Out [18]= (a - 2 n - 2 a x) F[x] + (-2 a + 2 x) F' [x] + F' ' [x] == 
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In [19]:= SimpleDE[Exp[a x~2] *Erf c [n,x] ,x] 

2 2 2 

Out [19]= (-2 a - 2 n - 4 a x + 4 a x ) F[x] + (2 x - 4 a x) F' [x] + 

> F" [x] == 

In [20] := SimpleDE [Abramowitz [n,x] ,x] 

(3) 

Out [20]= 2 F[x] + (1 - n) F" [x] + x F [x] == 
(see § (26.2.41)). 

We note that the algorithm obviously works for antiderivatives. An example of that type is 
Dawson's integral (see e. g. B (7.1.17)) for which we get the differential equation 

In [21] := SimpleDE [E~ (-x~2) *Integrate [E~ (t~2) ,{t,0,x}] ,x] 
Out [21]= 2 F[x] + 2 x F'[x] + F" [x] ==0 

For the Struve functions, our algorithm generates the differential equations 

(n 2 + n 3 + x 2 - nx 2 ) H n (x) + x (x 2 - n - n 2 ) W n (x) + (2 - n) x 2 H'^(x) + x 3 H^'(x) = , 

and 

(n 2 + n 3 — x 2 + nx 2 ) L ra (x) — x (x 2 + n + n 2 ) ~L' n (x) + (2 — n) x 2 L"(x) + x 3 L"'(x) = , 

that are the homogeneous counterparts of the differential equation (12.1.1) in Q|. 
Finally we give examples involving hypergeometric functions: 

In[22] := SimpleDE [Hypergeometric2Fl [a,b,c,x] ,x] 

Out [22]= a b F[x] + (-c + x + a x + b x) F' [x] + (-1 + x) x F" [x] ==0 

In [23] := SimpleDE [Hypergeometric2Fl [a,b,a+b+l/2,x] ~2,x] 

Out [23]= 8 a b (a + b) F [x] +2 (-a -2a -b-4ab-2b +x+3ax+ 

2 2 

> 2a x + 3bx + 8abx + 2b x)F'[x] + 

> 3x(-l-2a-2b + 2x + 2ax + 2bx)F"[x] + 

2 (3) 

> 2 (-1 + x) x F [x] == 



Here the last function considered 2F1 



- - 2 

a b 



x is the left hand side of Clausen's for- 



a+6+1/2 

mula ( |3l| ) that we will consider again in § |9[ 

Now we investigate the case that a derivative rule and a differential equation are given, and 
show that these two imply the existence of a recurrence equation: 
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Algorithm 4 If a family /„ is given by a derivative rule of type ([l]) and a differential equation 
of type (||), then it forms an admissible family for which a recurrence equation can be found 
algorithmically. 

Proof: We present an algorithm which generates a recurrence equation for f n : Iterative differen- 
tiation of the derivative rule (|l|) with the explicit use of (jlj) at each step yields 

M 

fn\x) = Y, r l( n > x ) fn-k{x) 

k=0 

with rational functions r° k . The substitution of these derivative representations in the differential 
equation gives the recurrence equation searched for. □ 

As an example we consider the Airy functions Ai n , again, for which we have the derivative rule 

© 

Ai n {x) = Ai n+ i (x) 

and the differential equation (|30| ) 

Ai"' (x) - x Ai' n (x) - (n + 1) Ai n (x) = . 

Differentiating the derivative rule successively and substituting the resulting expressions into the 
differential equation immediately yields the recurrence equation (|20|), again. 

If this family, however, is given by the backward derivative rule (compare (pC|)) 

Ai' n (x) = x Ai n _i (x) + (n - 1) Ai n _ 2 (x) , 

then differentiation yields 

Ai^ (x) = Ai n _i (x) + x A4_! (x) + (n - 1) Ai^_ 2 (?) 

= Ai n _i (x)+x(x Ai n _ 2 (x) + (n-2) Ai n _ 3 (x)J + (n-l)(x Ai n _ 3 (x) + (n-3) Ai n _ 4 (x)) 
= Ai n _i (x) + x 2 Ai n _ 2 (x) + (2n - 3) x Ai n _ 3 (x) + (n 2 - 4n + 3) Ai n _ 4 (x) . 

After a similar procedure we get 

Ai^' (x) = 3x Ai n _ 2 (x) + (x 2 + 3n - 5) Ai n _ 3 (x) + (3n - 6) x 2 Ai n _ 4 (x) 

+(3n 2 - 15n + 15) x Ai n _ 5 (x) + (n 3 - 9n 2 + 23n - 15) Ai n _ 6 (x) , 

and the substitution into the differentiation equation gives finally 

(n - 1) Ai n (x) - x 2 Ai n _i (x) + (4 - n) x Ai n _ 2 (x) + (3n - 5 + x 3 ) Ai n _ 3 (x) 
+(3n-6)x 2 Ai n _ 4 (x) + (3n 2 - 15n+15) x Ai n _ 5 (x) + (n 3 -9n 2 + 23n-15) Ai n _ 6 (x) = , 

a recurrence equation of order 6 rather than the minimal order three. This shows, that, in general, 
the order of the resulting recurrence equation is not best possible. 

Algebraically spoken, our result tells that if {fn \ j S No} has finite dimension, and if f n is 
an element of the linear space V spanned by a finite number of the functions {f n ±k}, then the 
space generated by all of {/ n ±fc} is of finite dimension, too. In contrast to Theorem |l], however, 
the dimension of this space generally may be higher than the dimension of V. This shows the 
advantage of the use of admissible families. 

As a further result of this section we note that using our general procedure developed in || we 
have 
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Algorithm 5 (Find a Laurent-Puiseux representation) Let / be a function that is built 
from the functions expx, lnx, sinx, cosx, arcsinx, arctanx, and any other functions that are 
embedded into admissible families, with the aid of the following procedures: differentiation, anti- 
differentiation, addition, multiplication, and the composition with rational functions and rational 
powers. 

If furthermore / turns out to be of rational, exp-like, or hypergeometric type (see then a 

oo 

closed form Laurent-Puiseux representation f(x) = a k % k can be obtained algorithmically. □ 

k=ko 



We remark that there is a decision procedure due to Petkovsek |19| to decide the hypergeometric 
type from the recurrence equation obtained. 

With Algorithm ||, it is possible to reproduce most of the results of the extensive bibliography 
on series [j?]], and to generate others. As an example we present the power series representation of 
the square of the Airy function: 

In [24] := PowerSeries [AiryAi [x] "2 ,x] 

lk k 1 + 3 k 
(-) 27 x (2k)! 
9 

Out [24]= Sum[-( ), {k, 0, Infinity}] + 

Sqrt[3] Pi k! (1 + 3 k) ! 

k 3 k 1 
12 x Pochhammer [- , k] 
6 

> Sum[ , {k, 0, Infinity}] + 

1/3 2 2 

3 3 (3k)! Gamma [-] 
3 

1/3 k 2 + 3 k 5 

2 3 12 (1 + k) x Pochhammer [- , k] 

6 

> Sum[ , {k, 0, Infinity}] 

1 2 

(3 + 3 k) ! Gamma [-] 
3 

Note that, moreover, this technique generates hypergeometric representations, whenever such rep- 
resentations exist. The above example, e. g., is recognized as the hypergeometric representation 



%2 1 ( V 6 

Al(x) 3 4 /3 T (2/3) 2 li?2 I 1/3 2/3 



^- & ( 1/2 

V3tt I 2/3 4/3 



9 y 3 2 / 3 T (1/3) 2 1 2 I 4/3 5/3 



9 , 



As soon as a hypergeometric representation is obtained, by Theorem ^ derivatives rules and re- 
currence equations with respect to all parameters involved may be obtained. As an example, we 
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consider the Laguerre polynomials: The power series representation for the Laguerre polynomial 
L { n\x) that our algorithm generates corresponds to the hypergeometric representation 



L { n a) (x) 



n + a 
n 



1F1 



—n 
a + 1 



from which by an application of Theorem |3] we obtain the derivative rule 

-n + 1 



dx n K ' \ n J x \ \ a + 1 



-(n + a) ( n - 1 + a 
n - 1 



1^1 



-(n-1) 
a + 1 



x + 



—n 
a + 1 

n ( n + a 



x 



x V n 



— n 
a + 1 



= ^(-( n + a)L { -:l l {x)+nL^{x)) , 
i. e. again, but we are also led to the derivative rule with respect a: 



n + a\a 
n )x[ lFl 



a n + a / n + a- 1 
x a V n 



—n 
a 



—n 
a + 1 



x 



lF! 



—n 
a 



x 



a r n+ a \ 
x V n ) 



— n 
a + 1 



= }-({n + a)L^- 1) {x)-aL^{x)) . 

A further application of Theorem [3| yields the recurrence equation 

1 + a ( 

F a+ i = ; — — - a F a _i + (a + x) F c 



(1 + a + n) x 



for F a := 



— n 
a + 1 



xj with respect to a, and the use of the algorithm for the product ( fllq 

Theorem 3 (d), J2^], p. 342, and J2l||, Maple function rec*rec), applied to -- 
generates dT9|), again. 



n + a 
n 



■F n 



9 Algorithmic verification of identities 

On the lines of |^7| we can now present an implementable algorithm to verify identities between 
expressions using the results of the last section. 

Algorithm 6 (Verification of identities) Assume two functions f n (x) and g n (x) are given, to 
which Algorithm || applies. Then the following procedure verifies whether /„ and g n are identical: 

(a) del :=SimpleDE(f ,x): 

Determine the simple differential equation del corresponding to f n . 

(b) de2:=SimpleDE(g,x): 

Determine the simple differential equation de2 corresponding to g n . 
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(c) (Different differential equation implies different function) If del and de2 have the 
same order, then 

- if they do not coincide besides common factors, i. e. have rational ratio, then f n and g n 
do not coincide; return this, and quit. 

- Otherwise f n and g n satisfy the same differential equation del of order I, say, and it 
remains to check I initial values. Continue with (e). 

(d) Let the orders of del and de2, i. e. 

I m 

J2pjA 3) =0 and X>^ fc) =0 
j=0 k=0 

(Pj (J = 0, . . ■ ,l),qk (k = 0, . . . , to) polynomials) are different, and assume without loss of 
generality that / > to. Then, differentiate de2 I — m times to get equations 

Sp--=Y,€9n k) =0 (p = m,...,l). 

k=0 

Check if there are nontrivial rational functions A p ^ (p = m, . . . ,1) such that a linear 
I 

combination S p is equivalent to the left hand side of del, i. e. is a rational multiple of 

p=m 

it. 

If this is not the case, then f n and g n do not satisfy a common simple differential equation, 
and therefore are not identical; return this, and quit. Otherwise they satisfy a common simple 
differential equation; continue with (e). 

(e) Let I be the order of the common simple differential equation for /„ and g n . For k = 0, ... ,1—1 
check if /i, (0) = g4 (0). (Note that by the holonomic structure the knowledge of the initial 
values (S) is sufficient to generate those.) These initial conditions may depend on n, and are 
proved by application of a discrete version of the same algorithm. If one of these equations 
is falsified, then the identity /„ = g n is disproved; return this, and quit. Otherwise, if all 
equations are verified, the identity f n = g n is proved. 

Proof: By a well-known result about differential equations of the type considered, the solution of 
an initial value problem 

I 

E Pk{x) fi k) (x) = , /«(0) = a k (k = 0, . . . , I - 1) 

k=0 

is unique. To prove that /„ and g n are identical, it therefore suffices to show that they satisfy a 
common differential equation, and the same initial values. This is done by our algorithm. □ 

For the example expressions 

/„(*) :=4~ 1/2) (*0 

and 

(-l) n 

9n(x) := ^2^ H 2n (v^) 
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we get the common differential equation 

2 n f{x) + (1 - 2x) f'(x) + 2xf"(x) = 0. 
Therefore to prove the identity 

4r 1/2) (*) = y> 2 n(v^) , 

(see e. g. ||], (22.5.38)), it is enough to verify the two initial equations / n (0) = g n (0) and /4(0) = 
g' n (0). To establish the first of these conditions, with Mathematica, e. g., we get 

In[25]:= eq = Limit [LaguerreL [n, -1/2 ,x] ,x->0] == 

Limit [ (-1 ) "n/ (n ! *2~ (2*n) ) *HermiteH [2*n , Sqrt [x] ] , x->0] 

1 

Pochhammer[l + n, -(-)] n 

2 (-1) Sqrt [Pi] 

Out [25]= == 

Sqrt [Pi] 1 

n! Gamma [- - n] 
2 

which is to be verified. In this situation, we establish the first order recurrence equations for both 
sides 

In[26] := FindRecursion [Limit [LaguerreL [n, -1/2 ,x] ,x->0] ,n] 
Out [26]= (-1 + 2 n) a[-l + n] - 2 n a[n] == 

In[27] := FindRecursion [Limit [(-1) ~n/(n ! *4~n) *HermiteH [2*n, Sqrt [x] ] ,x->0] ,n] 

Out [27]= (-1 + 2 n) a[-l + n] - 2 n a[n] == 

that coincide, so that it remains to prove the initial statement 

In [28] := eq / . n->0 

Out [28] = True 

and we are done. Similarly one may prove the second initial value statement /^(O) = g' n {0). 
Applying the same method, ( |T2| ) can be proved by the calculations 

In[29] := SimpleDE[(n+l)*Bateman[n+l,x]-(n-l)*Bateman[n-l,x] ,x] 

2 2 3 2 

Out [29]= (2n-2x + 4n x-4nx + x)F[x] + (-2nx + 2x)F'[x] + 

2 

> (2 n - x) x F" [x] == 
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In [30] := SimpleDE[2*x*D[Bateman[n,x] ,x] ,x] 

2 2 3 2 

Out [30]= (2n-2x + 4n x-4nx + x)F[x] + (-2nx + 2x)F'[x] + 



> (2 n - x) x F" [x] == 

and using the initial values F n (0) = and i^(0) = -2 (see fj], (11)). 
Also, one can prove Clausen's formula 



2-^1 



a b 
a+b+1/2 



x 



3F2 



2a 2b a + b 
a+b+1/2 2a+2b 



generating the common differential equation 



8a 6 {a + b) f(x) + 

2{-a-2a 2 -b - Aab-2b 2 + x + 3ax + 2a 2 x + 3bx + 8abx + 2b 2 x) f'(x) + 
3x{-l-2a-2b + 2x + 2ax + 2bx) f"{x) + 
2 (-1 + x) x 2 f"'(x) = 

for both sides of (|3l|), or other hypergeometric identities like the Kummer transformation 



1F1 



x\=e x xFx 



or like 



and 



)Fi ( a J x) ■ Fi ( b J x) = 2 F 3 ( 



b — a 
b 



a+b a+b— 1 
2 2 

a b a + b — 1 



Ax 



2^3 



a b — a 
b k b+i 

u 2 2 



X 

T 



\ b X )' lFl {i 

corresponding to the Kummer differential equation 

af(x)-(b-x)f'(x)-xf"(x) = 0, 

and to 

(l-a-b) (a + b) f(x) + {-ab + a 2 b + ab 2 -2x-Aax-Abx} f'{x) + 

+ (a + a 2 + b + 3ab + b 2 -Ax) xf"(x) + 2 (1 + a + b) x 2 f"'{x) + x 3 f""(x) = 0, 



and 



Aa (a-b) x f(x) + (b - 3 b 2 + 2 b 3 - x 2 - 2 bx 2 ) f'{x) 

+x (-6+5 6 2 -x 2 ) f"(x) + (l + Ab) x 2 f"'{x) + x 3 f""(x) = 0, 
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respectively. 

Note that one can also reverse the order of the algorithm, i. e. first find common recurrence 
equations for f n and g n with respect to n, and then check the initial conditions (depending on 
x) with the aid of differential equations. This method should be compared with recent results of 
Zeilberger 

Moreover the given algorithm is easily extended to the case of several variables, if the family 
given forms an admissible family with respect to all of its variables, i. e. for each variable exists 

- either a simple recurrence equation (corresponding to a "discrete" variable), 

- or a simple derivative rule (corresponding to a "continuous" variable), depending on shifts 
with respect to one of the discrete variables. 

Note, however, that (for the moment) the algorithm only works if / and g are "expressions", and 
no symbolic sums, derivatives of symbolic order, etc. occur. In the next sections, we will, however, 
extend the above algorithm to these situations. 



10 Algorithmic verification of Rodrigues type formulas 

Here we present an algorithm to verify identities of the Rodrigues type 

g(n,x) = /( n )(n, x) (f,g functions , n symbolic) . 

This algorithm, however, does only work if the function / is of the hypergeometric type. On the 
other hand, for most Rodrigues type formulas in the literature, see e. g. [||, this condition is valid. 
The procedure is based on the following 

Algorithm 7 (Find differential equation for derivatives of symbolic order) Let / be of 

the hypergeometric type, i. e. there is a Laurent-Puiseux type representation f(n,x) = J2k a k xk - 
Then there is a simple differential equation for g(n,x) := f^ n \n,x) which can be obtained by the 
following algorithm: 

(a) del :=SimpleDE(f ,x): 

Calculate the simple differential equation del of /, see Algorithm^ 

(b) rel :=DEtoRE(del ,f ,x,a,k): 

Transfer the differential equation del into the corresponding recurrence equation rel for a k , 
see [01, §6. 



(c) If rel is not of the hypergeometric type (or is not equivalent to the hypergeometric type |19| ), 
then quit. 

(d) re2 : =SymbolicDerivativeRE(rel , a,k,n) : 

Otherwise set ct '■= {k + l) n a,k+ n - Bring rel into the form 

afc+m = R{k) a k , 

rational R, and calculate the hypergeometric type recurrence equation re2 

(k + n + l) m 

Ck+m = /; , 1N R{k + n) c k (32) 

{k + l) m 

for c fc . 
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(e) de2:=REtoDE(re2,a,k,G,x): 

Transfer the recurrence equation re2 into the corresponding differential equation de2 for the 
nth derivative g(n,x) := f <yn \x) of /, see |9j, § 11. 

Proof: Parts (a), (b) and (e) of the algorithm are described precisely in |]]. Now, assume, 
g(n,x) = f( n \n,x), and that / has the representation f(n,x) = J2 k a k x k . Then we get 

Y^c k x k = g(n,x) = f {n) (n,x) = ^0 + 1 ~ n) n a k x k - n = ^(k + l) n a k+n x k . 

Therefore we have c k = (k + l) n a k+n , and we get the recurrence equation 

c k+m = (k + m + l) n a k+n+m = (k + m+ l) n R(k + n) a k+n 

{k + m + l) n + m + n)! k\ 

= ruj-u R(k + n)c k = R{k + n)c k 

{k + l) n (k + m)\ (k + n)\ 

{k + m + n)\ k\ . . (k + n + l) m 

= — ( 1 I \ — 77— — r-.R{k + n)c k = . R(k + n)c k , 

(k + n) (k + m)\ (k + l) m 

and hence (^), for c k . This finishes the proof. □ 
As a first example we consider the identity 

erfCn(x) = 2 "n! V 6 " erfCX ) 

(see e. g. Q, (7.2.9)), or equivalently 

(-1)™ 2 n n\ e x2 erfc n (x) = (e* 2 erfcx) . (33) 

Algorithm yields step by step 

In[31] := del=SimpleDE[E-(x~2)*Erfc[x] ,x] 

Out [31]= -2 F[x] - 2 x F'[x] + F" [x] ==0 

In[32]:= rel=DEtoRE[del,F,x,a,k] 

0ut[32]= -2 (1 + k) a[k] + (1 + k) (2 + k) a[2 + k] ==0 

In[33] := re2=SymbolicDerivativeRE[rel,a,k,n] 

2 

Out [33]= -2 (1 + k + n) a[k] + (2 + 3 k + k ) a [2 + k] ==0 
In[34]:= de2=REtoDE[re2,a,k,G,x] 

Out [34]= -2 (1 + n) G[x] - 2 x G' [x] + G" [x] == 
thus finally the differential equation 

-2{l + n)g(x) -2xg{x) +g"(x) =0 
for the function -jj^ (e x2 erfc x^j , which also can be obtained by the single statement 
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In [35] := RodriguesDE [E~ (x~2) *Erfc [x] ,x,n] 

Out [35]= -2 (1 + n) F[x] - 2 x F' [x] + F" [x] == 

For the left hand term of (|3^) we get 

In [36] := de3=SimpleDE [E~ (x~2) *Erf c [n,x] ,x] 

Out [36]= -2 (1 + n) F[x] - 2 x F' [x] + F" [x] == 

i. e. the same differential equation. 

As next example we consider the Rodrigues type id 
rewrite it as 



Our implementation yields 

In [37] := RodriguesDE [E~ (-2x) *x~ (n-l) ,x,n] 

Out [37]= 2 (1 + n) F[x] + 2 (1 + x) F' [x] + x F" [x] == 

In [38] := SimpleDE [E~ (-x) /x*Bateman [n,x] ,x] 

Out [38]= 2 (1 + n) F[x] + 2 (1 + x) F' [x] + x F"[x] == 

Algorithm shows the applicability of Algorithm |6| if in the expressions involved Rodrigues type 
expressions occur, as soon as we can handle the initial values. Since in Algorithm |7] the function 
/ is assumed to be of hypergeometric type, this, however, can be done by a series representation 
using Algorithm || if / moreover is analytic, and if the function / of Algorithm ^ does not depend 
on n: In this case Algorithm [5] generates the generic coefficient a/% of the series representation 



f{x) = J2 Q-kX k , and therefore we get the initial values by Taylor's theorem: 
fc=o 





(34) 



oo 




In our first example we conclude 



In [39] : 



PowerSeries [E~ (x~2) *Erf c [x] ,x] 



2 k 



x 



Out [39] 



Sum[ , {k, 0, Infinity}] + 

k! 



k 1 + 2 k 



-2 4 x 



k! 



> 



Sum[ 



, {k, 0, Infinity}] 



Sqrt[Pi] (1 + 2 k) ! 
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so that the first initial condition for identity (|33|) is given by the calculation (see pf (7.2.7)) 



(2*0 



^J-_^ = (_l)»2»n! erfc n (0) = ^ (e* erfcx) (0) J 



if n = 2fc (A; G N ) 



if n = 2fc+l (fe GN ) 



and the second one is established similarly. 

To identify the first initial values of our second example, we proceed as follows: The left hand 
side of (H) yields 

Jim '±^— F 'Jx) = n! lim = n\F'(0) = -2n\ (35) 

x->0 x x->-0 X 



, /^x \ n\ ii k = n 
otherwise 



(see |14| ) (11)), whereas from the identity 

(x") W (0) 

and Leibniz's formula we derive for the right hand side 



\k=0 



( n ^)(n-l)! (e- 2 j'(0) = -2n! 



in agreement with 

It is easily seen that we can always identify the initial values algorithmically by the method 
given if /(n, x) = w(x) X(x) n with a polynomial X, i. e. is of the form (Q). 
These results are summarized by 

Algorithm 8 (Verification of identities) With Algorithms |6| and fj] identities involving Ro- 
drigues type expressions can be verified if only symbolic derivatives of hypergeometric type 
analytic expressions / occur that have the form f(n,x) = w(x) X{x) n for some polynomial X. 

11 Algorithmic verification of formulas involving symbolic sums 

In this section we study, how identities involving symbolic sums can be established. The results 
depend on the following algorithm (compare [pi]], Maple function cauchyproduct): 

Algorithm 9 (Find recurrence equation for symbolic sums) Let f n (x) form an admissible 

n 

family, and let s n (x) denote the symbolic sum s n (x) := J2 fk(x)- Then the following algorithm 

k=0 

generates a recurrence equation for s n : 

(a) re : =FindRecursion(f ,k) : 

Calculate the simple recurrence equation re of see Q, §11. 

(b) del :=REtoDE(rel ,f ,k,F,z): 

Transfer the recurrence equation re into the corresponding differential equation del valid for 

oo 

the generating function F(z) := J2 fk( x ) z k , see ||, §11. 

k=0 
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(c) de2:=F(z)+(z-l)*F> (z)=0: 

Let de2 be the differential equation corresponding to the function 

OO OO -I 

G(z) :=E^ fe = E^ = T37- 

k=0 k=0 

(d) de:=ProductDE(del,de2,F,z): 

Calculate the simple differential equation de corresponding to the product H{z) := F{z) G(z), 
see [jl5f| , Theorem 3 (d). This differential equation has the order of del. 

(e) re:=DEtoRE(de,F,z,s,n): 

Transfer the differential equation de into the corresponding recurrence equation re for the 
coefficient s n of H(z), see ||, § 6. 

Proof: Parts (a), (b) and (e) of the algorithm are described precisely in The rest follows from 
the Cauchy product representation 

oo/n \ oo / n \ 

n=0 Vfc=0 / n=0 \fc=0 / 

of the product function F(z) G(z). □ 

n i 

As an example we consider the sum J2 L>k (x). We get stepwise: 

fc=o 

In[40] := re=FindRecursion[LaguerreL[k, alpha, x] ,k] 

Out [40]= (-1 + alpha + k) a[-2 + k] + (1 - alpha - 2 k + x) a[-l + k] + 

> k a[k] == 

In[41] := del=REtoDE[re,a,k,F,z] 

2 

0ut[41]= (-1 - alpha + x + z + alpha z) F[z] + (-1 + z) F' [z] == 
In[42] := de2=F [z] + (z-l) *F' [z]==0; 
In [43] := de=ProductDE[del,de2,F,z] 

2 

Out [43]= (-2 - alpha + x + 2 z + alpha z) F[z] + (-1 + z) F' [z] == 
In[44] := DEt oRE [de , F , z , s , n] 

Out [44]= (2 + alpha + n) s [n] + (-4 - alpha - 2 n + x) s[l + n] + 

> (2 + n) s[2 + n] == 

or by a single statement 
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In[45] := re=SymbolicSumR.E[LaguerreL[k, alpha, x] ,k,n] 

Out [45]= (2 + alpha + n) a[n] + (-4 - alpha - 2 n + x) a[l + n] + 

> (2 + n) a [2 + n] == 

and substituting n by n — 2 

In [46]:= Simplify [re /. n->n-2] 

Out [46]= (alpha + n) a [-2 + n] + (-alpha - 2 n + x) a[-l + n] + n a[n] ==0 

On the other hand, the calculation 

In [47] := F indite cur si on [LaguerreL [n, alpha+1 ,x] ,n] 

Out [47]= (alpha + n) a [-2 + n] + (-alpha - 2 n + x) a[-l + n] + n a[n] ==0 

shows that the left and right hand sides of the identity 

n 

£4 a) (*)=4 a+1) (z) (36) 

k=0 



(see e. g. [24], VI (1.16)) satisfy the same recurrence equation. 

In our example identity two initial values remain to be considered 

Z>) ( x ) = Z> +1 ) ( x ) = l and 4 a) (x) + L<f ] (x) = L^ +1) (x) = 2 + a - x 

that trivially are established. 

Thus Algorithm || shows the applicability of Algorithm |6| if in the expressions involved symbolic 
sums occur. This is summarized by 

Algorithm 10 (Verification of identities) With Algorithms || and || identities involving sym- 
bolic sums can be verified. □ 

We like to mention that the function FindRecursion is successful for composite f n as long as 
recurrence equations exist and are applied recursively. Here obviously no derivative rules are 
needed. 

We note further that as a byproduct this algorithm in an obvious way can be generalized to sums 

n 

a-k b n -k of the Cauchy product type. As an example, the algorithm generates the recurrence 

fc=o 

equation 

2(1 + 2n) s n - (1 + n) s n+1 = (37) 
i l 



for Sn := E (I) =nl 2 E^ ¥ 



x ,. , ^MF' 

k=0 n ' k=0 
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In [48] := re=ConvolutionRESum[l/k! ~2, 1/k! ~2,k,n] 



3 

Out [48]= 2 (1 + 2 n) a[n] - (1 + n) a[l + n] == 



In[49] := ProductRE[re,FindRecursion[n! ~2,n] ,a,n] 



Out [49]= 2 (1 + 2 n) a[n] + (-1 - n) a[l + n] == 



compare 



Algorithm || may further be used to find a closed form representation of a symbolic sum in case 
the resulting term is hypergeometric: 



n 



Algorithm 11 (Closed forms of hypergeometric symbolic sums) Let s n := fk be a 



hypergeometric term, i. e. ^y^- be a rational function, then the following procedure generates a 
closed form representation for s n : 

(a) re : =SymbolicSumRE(f ,k,n) : 

Calculate the simple recurrence equation re of s n using Algorithm ||. 

(b) If re is of the hypergeometric type, then solve it by the hypergeometric coefficient formula, 
else apply Petkovsek's algorithm to find the hypergeometric solution s n of re. □ 

This result should be compared with the Gosper algorithm O. Our procedure is an alternative 
decision procedure for the same purpose. Note that from the hypergeometricity of s n the hyper- 
geometricity of fk follows ||, so that the first step of Algorithm ^ leads to a simple first order 
recurrence equation. 

Applying our algorithm to our example case s n = [ u ) > we § e t from ([37]), and the initial 



fc=o 




value so = 1 the representation 




On the other hand, for s n = J2 ( u ) > our procedure gives 

fc=o v fc ' 

In [50] := re=ConvolutionRESum[l/k! "3, 1/k! ~3,k,n] 




2 

Out [50]= 8 a[n] + (1 + n) (16 + 21 n + 7 n ) a[l + n] - 



> 



5 

(1 + n) (2 + n) a [2 + n] == 



In[51] := ProductRE[re,FindRecursion[n! ~3,n] ,a,n] 



2 2 2 

0ut[51]= 8 (1 + n) a[n] + (16 + 21 n + 7 n ) a[l + n] - (2 + n) a[2 + n] == 



> 
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and by Petkovsek's algorithm it turns out that s n is no hypergeometric term. 
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